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RESEARCH MEMORANDUM 


THEORETICAL STUD! OF THE TUNNEL-BOUNDARI LIFT 
INTERFERENCE DUE TO SIiOTTED WALLS IN THE 
PRESENCE OF THE TRAILING-VORTEX 
SYSTEM OF A LIFTING MODEL 
By Clarence W. Matthews 

SUMMARY 


The equations presented in this paper give the interference on the 
trailing -vortex system of a uniformly loaded finite- span wing in a 
circular tunnel containing partly open and partly closed walls, with 
special reference to symmetrical, arrangements of the open and closed 
portions. Methods are given for extending the equations to include 
tunnel shapes other than circular. The rectangular tunnel is used to 
demonstrate these methods. The equations are also extended to non- 
uniformly loaded wings. 

An analysis of the equations for certain configurations has shown 
that; (l) only a small percentage of slot opening is required to give 
zero interference conditions if the tunnel contains four or more slots j 
( 2 ) in the configurations studied, the ratio "between the slotted-tunnel 
interference and the closed-tunnel interference at the center of the 
tunnel is approximately constant for various model spans; and ( 3 ) timnels 
containing an odd number of slots or nonsynanetrical slot arrangements 
cause an additional rolling moment or a cross flow on the wings, or both. 


INTRCDUCTION 


In a study of solid -blockage Interference (see ref. 1), it has been 
shown that tunnels containing mixed boundaries, that is, partly open and 
partly closed walls, will eliminate or greatly reduce such interference. 
Such tunnels have been shown to be very useful test equipment (ref. 2). 
Since the slotted tunnel coixC igviration required to eliminate solid 
blockage may not eliminate lift interference, it is necessary to study 
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the Interference on the trailing vortices of a lifting model in order 
to iuftlce the necessary corrections to the lift characteristics of a model. 

The problem of one or two slots tias been treated by varioxzs authors 
(see refs, 3 to 6). The case of more than two slots has also been treated 
in references ^ and 8 . Reference 7 treats only small wings in circular 
tunnels, and reference 8 treats only the case for large numbers of evenly 
spaced slots. 

The purpose of this report is to present equations which express 
the tunnel- wall interference due to mixed open and closed boundaries in 
the presence of the trailing- vortex Systran of a finite-span lifting model 
at subsonic velocities. Special attention is given to test sections in 
which ttie slots ane symmetrically located with respect to both axes. 

Various extensions of the theory have been made, following in a 
general fashion the methods of reference 7 * These extensions include 
the effects of wing span, slot configurations, interference at points 
near the center, nommiform loading, and methods for calculating the 
interference in tunnels of other thm clrculer cross section. 

Numerical calculations of the interference characteristics of several 
symmetrical cases are presented and are used to show the properties of 
the interference of circular tunnels containing 1, 2, k, 8, and 12 slots 
symmetrically located with respect to the x- and y-axes, and a square 
tunnel containing 8 slots symmetrically located in the top and bottom 
walls. 


SYMBOLS 






ctn^ Pn 

b 

C 

Cl 



>. constants 

half- span of wing 

txmnel cross-sectional area 

lift coefficient of model 


c 
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Cd 

ACd 

ds" 

S 

k 


I 

m 

M 


n,s,q. 


P 

r 


R.P. 


I.P. 


S 

u 

V 

V 

x,y 

z 

dw 

dz 

r 

7 


drag of model 

drag increment due to interference 
tangential line element 
multiple -product index 

quality factor, the ratio at any given point of the inter 
ference of a slotted tunnel to the interference of a 
closed tunnel with the sane cross section 

half -span length of wing in ^-plane 

number of slots in tunnel 

Mach number 

summation indices 

special number defined in equation ( 2 ) 

radial distance of point from coordinate center 

real part of canplex function 

imaginary part of complex fvinctlon 

model or wing area 

X- component of velocity 

y-con^onent of velocity 

velocity vector 

coordinates of rectangular system (fig. 1) 
coti^lex variable, x + iy 

coD^lex velocity, u - iv 

circulation about a point, positive in counterclockwise 
direction 

function of Sji and b defined in equation (32) 


k 
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8 

a 

5 

e 

■Q 

\ 

9 

^cp 

dn 


proportionality factor 

Interference correction angle to be ewided to measured 
angle of attack 

strength of a soiirce 

congjlex variable in transformed plane 

angular coordinate of polar system (fig. 1 ) 

complex function of z (z = e^'®) 

parameter defined in eq.uatlon (105) (also see par. 14.8^ 
ref. 10 ) 

pertizrbation potential 

derivative of potential in direction normal to a given line 



THEORT OF LIFT-VOHEEX 


nr SLOTTETD TIMMELS 




General Analysis 


• Theoretical boundary conditions of flow about trailing vortices in 
a tunnel containing mixed open and closed portions .- The equations for 
the interference on the lift of a model due to mixed open and closed 
tunnel walls in the presence of a traillng-vortex system may be obtained 
by considering the same two-dimensional approximation of the flow field 
that is used in reference 9 * The conditions of this two-dimensional 
approximation may be briefly stated as; (l) the tunnel and its boundaries 
extend from a point an infinite distance upstream of the model to a point 
an infinite distance downstream of the model; (2) the velocities Induced 
by the trailing vortices in the cross-sectional plane located at the model 
are one-half of those induced in a far-downstream cross-sectional plane; 
( 3 ) the induced- velocity flow field in the far- downstream section may be 
treated with two-dimensional methods; (4) the boundary condition which 
must be satisfied at a solid portion of a tunnel wall is that the flow 
must be tangential to it^ or 
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(5) the condition which must he satisfied at an open portion of a tvinnel 
wall is that the potential over that portion must he constant or the flow 
must he normal to itj (6) no slngular’ities other than the trailing vortices 
can exist within the hotmdaries of the tunnel; and ( 7 ) the constant 
potential in every slot must he equal to that in every other slot. The 
final condition is required because the pressure is the same at every 
slot and may he shown hy the following considerations. Since the pres- 
sure is constant over the entire region outside the tunnel there will 
he no pressure differences between the slots at a point far upstream^ 
and hence there can he no flow between the slots due to external influ- 
ences. Also, since this point is too far from the model to he influ- 
enced hy it, there can he no flow due to the model. As there is no flow 
between the slots, no potential gradient can exist between the slots; 
therefore, the potentials in all the slots must he the same. It there- 
fore follows hy condition (5) that in the far-downstream position, the 
constant potential in eveiy slot is equal to that in every other slot. 

Velocity fields in circular slotted tunnels .- The previoiisly stated 
boundary conditions may he satisfied hy using complex velocity functions 
rather than complex potential functions. This is done hy selecting a 
complex velocity fvmetion that has singularities at the wing tips and 
has a flow direction either nonnal or tangential to the tunnel walls. 

If this function is multiplied hy another complex f^mctlan whose value 
on the tunnel wall changes firom all real to all imaginary, or the 
opposite, at each slot edge, then the flow will he rotated 90 ° at those 
points so that the final flow of the prod;ict of the two functions will 
he normal to the watll on selected portions and will be tangential to 
the wall on the remaining portions. It may, however, he expected that 
the second function will Introduce singularities within the tunnel which 
are not permitted according to the stated boundary conditions, so that a 
third function which contains all the forbidden singularities must he 
\ised in such a fashion that it will cancel the forbidden singularities 
of the second function. 

In order to make up the first function, suppose that the two singu- 
larities at the w ing tips with their reciprocal singularities are written 

as l/(z2 _ h2)(l - z^h^), where z is the complex variable x + iy 
and h is the semlspan of the lifting wing. If this fxmctlon is examined 
hy letting z = e^®, where d is a polar angle (fig. 1), it will he found 
that the flow may he made normal to the walls if the factor z is 
Included in the nimnerator . This flow may also he made tangential to the 
wall hy multiplying hy the factor i. Thus, the first function may he 
written 


pz 

(z2 - h^j^l - z^-faS) 


( 2 ) 
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The symbol p may be chosen to be either 1 or 1, depending on whether 
normal or tangential flow is required at the walls. 


The second function may be developed by considering the square root 
of a frinction which is real on the wa.11. and changes sign at each slot edge 
so that the square root of the function changes frcaa real to imaginary at 


each slot edge. Such a function may be escpressed by ycot - cot ^ 
where z = e^'® and iS is in general complex. Examination of this func- 


tion shows that it becomes 


1 9 

cot - cot 


on the timnel wall, where 


9g is the polar angle of a slot edge (see fig. 1). The term under the 


radical changes frcm positive to negative at 9 = 9g so that the fimc- 
tion changes from re'al to imaginary at the point 0g on the tunnel wall. 
If a number of these functions having different 6g's marking the transi- 
tion from open to closed sections of the tunnel wall are multiplied 
together, it may be seen that the product will change sign at each value 
of ©g so that the fimction will be real on alternate sections. This 

use of cot 5- also sviggests the use of other ccorplex trigonometric func- 

tions such as cos d, sin d, and tan The sine and cosine terms will, 

on examination, be fotond to introduce two slot edges for each value of 0g, 
rather than one. These two slot edges will be found to be located symmet- 
rically with respect to the x-axls for the cosine term and to the y-axls 
for the sine term. Thus, the use of these functions is siiggested for 
symmetrical slot configurations. 


Since any of these fimctions inay Introduce singularities within the 
tunnel, it is necessary to examine them for such singularities . The four 
functions may be written as the following multiple products. 


<1 

n \jcOB ■5 - cos 0g s ^COS ■d - cos ©2_ \jc03 ^ - cos ©2 
g=l 


. ycos d - cos 0g 


(3) 



^sln d - sin 0g 


ih) 


q. 


n 

g=i 




( 5 ) 
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and 

n ^tan J - tan ^ (6) 

When z = is substituted for ■6, the above equations become, 

respectively. 


2 < 1/2 + 1 - 2z cos 0^ 

g=l 


S 


(7) 


q. 


^-q/2^-q/Z^-q/2 1 - 2iz sin 6g 


( 8 ) 


(z - /i(z + 1) - (z - l)cot ^ 

g=l 


(9) 


(z + J{z - 1) - i(z + l)tan (lO) 

g=l '' 

Examination of the f\mctions ( 7 ) to (lO) shovs that they contain 
forbidden singularities at the following valties of z: the cosine and 

sine terms at z = 0, the cotangent term, at z = 1, and the tangent 
term at z = -1. Several functions which contain singularities identical 
to those appearing in f motions ( 7 ) to (lO) and which are resQ. on the 
t unn el are 



(an cos n^a + bn sin na) 


( 11 ) 
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which has a singularity at z = 0; 


q/2 

^ (on + i|3n)cot“^|) 


( 12 ) 


which has a singularity at z = Ij and 


n=0 


(on + iPn)tan“^|^ 


(13) 


which has a singularity at z = -1; where an, an^ "bn^ Pn 

real constants. It may he seen now that hy dividing the trigonometric 
series hy the multiple product \rtiich contains the same singularities, 
the forbidden singularities will he canceled out of the final equation, 
leaving an eq'uation which has only the desired singularities within the 
tunnel. 


The reason that the multiple product is chosen for the denominator 
may he seen hy examining the flow about the slot edges. Since the flow 
must turn a sharp comer as it goes around the edge of each slot, the 
velocity must he infinite at that point. Since the multiple -product 
function becomes equal to zero at each slot edge, it must he placed in 
the denominator of the final can 5 )lex velocity function to insure the 
required infinite velocities. 


Before writing the final complex velocity function, it is observed 
that function (2) has a zero at the point z = 0 v^hich must also he 
removed hy either including an extra term in summation (U) or hy 
multiplying summation (l2) or (13 ) by a;j_ cos i8 + h^ sin d. With this 
note in mind, the final complex velocity functions may he written; 


dw 

dz 


pz ^ (an cos nd + hn sin nd) 
n®0 

q. 

_[ ^ ^COS d - cos 
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dw 

3^ = 




pz y (Sq cob ri'A + "bn sin nd) 
n=0 

q. ^ 

(z^ - b^)(l - z^b^) J_ _j_ Jain -9 - sin 0 g 

g=l 


(15) 


dw 
^ “ 


q/Z 


^ / \ 

pz(tt 2 ^ cos i 8 + bi sin d) j_ J_ (ocq + lgn)cot°( 
E^O 

f z^ - b^ ) fl - z^^) 1 r /cob — - cot 


(16) 


n n 

and an equation similar to (l 6 ), using tan ^ rather than cot 

In the application of these equations, it may be observed that each 
value of 0g introduces two slot edges into each of equations ( 1 I 4 -) 
and (l5) and only a single slot edge into eq.uation ( 16 ). Since each 
panel has two edges, q. must be equal to the number of panels m in 
equations (ll)-) (l5) and to twice the number of panels, or 2m, in 

equation ( 16 ) . Also, if equation (l^J-) or (15) is applied to tunnels 
containing an odd mmber of slots, it will be fo\md that the singularity 
arising from the multiple product will contain a term of the order I/ 2 , 
which cannot be removed by the summation. ]n order to remove this sin^- 
larlty, the summation of equation (ill-) or (15) must be rewritten as 


m-^l 

^ (a„ cos (17) 

This series will, upon examination, be found to contain a singularity of 
the order 1 / 2 , which will cancel the one-half -power term due to the exten- 
sion of the multiple product over an odd number of slots. 

In order to make the final application of equations (l^), (15), 
or ( 16 ) to a wind tunnel, the arbitrary constants of the simmiations must 
be evaluated. This evaluation can be effected by iising the remaining 
boundary conditions, which are that the potent ieils in each slot must be 
equal to each oldier and that only a vortex flow can exist about the 
singularities within the ttmnel. 
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The potential condition may he evalviated by the following line 
integral; 


p(X2 72) _ 

L’n) ° 

in which the limits of Integration terminate in the slot. One path which 
may be used is the streamline which flows along the panel. When this 
streamline is used as the path of integration, equation (l8) becomes 


''®n+l 

».®n 


V d0 = - cpa = 0 


(19) 


as the potential is the same over all slots. Since the velocity V will 
in general be complex, two equations exist for each panel. 

The strengths of the circulation about a pole and, of the soiirce at 
a pole are fixed by the relation. 


r + ia = ^ dz (20) 

where F is the circulation due to lift and a is the source strength 
and is equal to zero. Since equation (20) must be evaluated about each 
pole, it gives four equations which may be used to evaluate the constants. 

The set of boimdary equations (l9) and (20) may be considered as 
a set of simultaneous equations in the \mknown constants found in equa- 
tion (l4), (l5), or (l6) and is used to determine the values of these 
constants. If these constants are to be uniquely determined, there must 
be as many equations as there are consteints. An examination of equa- 
tions (19) and (20) shows that equation (. 19 ) gives 2m equations and 
that equation (20) gives four equations, making a total of 2m + 4 
equations which may be used to determine the \mknown constants. The 
number of constants which must be satisfied is determined frcm an exami- 
nation of equation (l6) . It is shown in the paragraph following equa- 
tion (l6) that q must equal 2m, so that the number of constants must 
equal 2m + 4; hence the number of equations and the number of ccaastants 
are equal, so that all the constants may be imlqviely determined, thus 
giving the problem a unique solution. 
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Derivation of Equations 

Slots symmetrically located, vlth respect to the x- and y-axes .- 
In the case of symmetry about the x- and y-axes, considerable simplifica- 
tion results. The slnpler equation (llf) may be Tised In place of equa- 
tion (l6). In equation (l 4 ), it may be observed that each valtie of ©g 
In the multiple products can be used to produce four changes of sign, 
provided these products are written as 

m/2 

n 

g=i 


m/2 

n 

g=i 

These changes of sign are observed at 0g, -6g, n - 0 g, and n + 0g, 

or one in each quadrant. Hence, if the values of 0g are known in one 
of the quadrants, the multiple product need be extended only over that 
quadrant, as al l the other slots will be automatically introduced. If 
the slots are evenly spaced and of equal width, function (21) or (22) 
can be replaced by 

^os^ 1= ^ - cos^ s (23a) 

or 

^sin^ ^ ■a - sin^ 5 0;j_ (23b) 

That (23a) or (23b) produces slots in the desired location may be seen 
by substituting 


^cos^’8 - cos^0| 


g 


/sin^’a - sin^0g 


(21) 


(22) 


i ,p = 0 

m ^ 


on the tunnel wall and 


0n = — - 

1 m. 


9 l 


where m is the n^miber of slots in the tunnel and n is any integer. 
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TaMng this modification into consideration, equation (ll)-) may he 
written with the suhstltutlon of z = e^'® as 



which may he reduced to 


dw 

dz 


f - 1 

2^ p 


n=>0 I_ 


m - 
« - n+1 , 


anz‘ 


,2n 


+ 1 


) - 


ibnz 


|-n.l 


(z^ - l) 


m/2 

(z^ - h^) (l - z^h^) [ _[ \jz^ - 2z^cos 20g + 1 

g=l 


(25) 


where the factor p is eqxial to 1 if slots intersect the x-axis and is 
equal to 1 if panels intersect the x-axis. If the slots are evenly 
spaced, and of equal width, the function 


/z2m - 2z“co8 mSi + 1 

may he substituted for the multiple product in equation (24) or (25). 

The constants of equation (25) must he evaluated hy using the bound- 
ary condition for the equality of the potential in each slot, equa- 
tion (19) ^ and for the circulation and nonexistence of sovirces, equa- 
tion (20). Equation (20) is evaluated hy applying Cauchy's integral 
theorem to equation (25) to obtain the required line integral of dw/dz, 
or 


2itip2 


n, - + 1 

2 — 


r + ia = 


n=i 


aj^h 2 + 1) - ihjib2 " - l) 


m 

j-n+l, 


m/2 


(26) 


2h(l - h^) J_ 2. l/h^ - 2h2cos 20g + 1 

g=l 


for the value of the Uae Integral about the positive pole b, and 




for the value of the line Integral about the negative pole -b. The circulation (eq,. 20) la 
determined by setting F equal to the real parts of equations (26) and (27) while the non- 
existence of sources is assured by setting the imaginary parts of equations (26) and (27) equal 

to zero. The condition on the potential, equation (19); is satisfied by letting z = e^® in 
equation (l^), thus obtaining the velocity along the panel streamline, which is a convenient 
streamline to use for the evaluation of equation (l9). Thus, eqmtlon (19) beccanes 



where the symbol 


indicates integration over a panel. Separate equatlcxiB are then 


obtained, erne for each panel. The constants and b^^ may now be evaluated by considering 
equations (26), (27)^ (28) as a set of simultaneous equations In and bj^. 


H 

u» 


NACA EM I53A26 



KACA EM L53A26 




In solving equations (26) to (26) for an. and tn, it is observed 
that for each a^; 


&n = On 



(29) 


with a corresponding relation between hn. and a new constant Pn.* With 

the substitution of these new constants On and Pn., eq.uation (25) may 
be written 


m/2 


dw 

dz 


pbr(i - b^) n - 2b^cos 20g + 1 r 


g=l 


m/2 


n(z2 - b2)(l - zV) - 2z2cos 20g + 1 

g=l 


o^z 

n=0 •— 


|-n+l/ V 

2 [z^°- + ij 




S-n+1 


(z2n 




(30) 


and eq,uations (26), (27), and (28) become, respectively. 




i 

P 


5 L- 


(b^^ + 1 


S-n+1 
2 


) - ip^-b (b^^ - 1) 


(31a) 






m 


±^(-•6)2 (31b) 


n=0 
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0 = 


m T 
2 *^ 


y (ot)2 C08 n© - iPn sin n©)d0 
n=0 

m/2 


■'panel 


(l - 2b^cos 20 + /cos® 0 - cos^0j 

g=l 


S 


( 31 c) 


vhere and 3j2 are the new constants. 

It can he seen hy examining equations (3I) that a number of the 
and equal to zero. This is shown hy first considering 

the nature of the simultaneous equations ( 31 )* The equating of the real 
and imaginary parts divides the entire set into two separate sets of 
simultaneous equations^ one of which involves o^'s only and the other 
3 j3^’s only. When p is chosen, one of the sets will he homogeneous 
whereas the other set will contain two equations with constant terms 
equal to -1. Thus^ one set of constants hecomes equal to zero whereas 
the other set will have values which may he dif ferment from zero. This 
remaining set of equations can he reduced hy considering the values of 
the integrals across any one of the panels and its counterpart which is 
symmetrically located with respect to the y-axis. It can he shown hy 
substituting ir - cp for 0 in these Integrals that they w in have the 
same absolute values and the same sign if n is even hut opposite signs 
if n is odd. The same featiire is also observed in ccKparlng equa- 
tions (31a) and (3lh); that is, if — _ n + 1 is odd, then the terms 

which, contain h to that power will have opposite signs whereas the 
remaining terms will have the same signs. It can now be seen that, hy 
adding or subtracting equations (31a) and (31b) and the equations for 
the integrals across each of the symmetrically located panels, it is 
possible to reduce the set of equations in c% (or in Pn) to two new 
sets, one of which contains the odd values of n and the other the even 
values of n. Again one of these sets is homogeneous, and one equation 
of the other set has a constant term equal to -1 so that a.1 1 the constants 
associated with either the odd values of n or the even values of n 
will be equal to zero and the other set will have values which can be 
different from zero. 

This analysis shows that one of four sets of constants <X2s, “gs+l^ 
P2s+1 occur, thus siaggesting fotn- different symmetrical 

slot configurations. The possible synmetrlcal configurations which can. 
occ\ir are (l) panels Intersecting both axes, (ll) slots intersecting both 
axes, (ill) a panel intersecting the x-axls and a slot intersecting the 
y-axis, and (IV) a slot Intersecting the x-axis and a panel intersecting 
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the y-axis. These conditions are found to occvir when the tunnel has 
i^-s slots with p equal to i (case I) or 1 (case II) or 2 ( 2 s + 1) 
slots with p equal to i (case III) or 1 (case IV), Since there are 
four sets of constants and four different symmetrical slot configura- 
tions^ it is to he expected that each set of constants can he associated 
with one of the symmetrical slot configurations. VThen equations (31) 
are set up for each of the symmetrical slot configurations^ it is found 
that the following associations exist. 

Case I, panels intersecting hoth axes: the oc^gs+l) values 

different from zero. 

Case II, slots intersecting hoth axes: the P( 2 s+l) values 

different frcxa zero. 

Case III, panel intersecting the x-axis and slot intersecting i?he 
y-axis: the set has values different from zero. 

Case IV, slot intersecting the x-axis and panel intersecting the 
y-axis : the hs set has values different from zero. 

It is now possible to write equations (30) and (3I) hy using only 
the constants which may have values other than zero. However, equa- 
tions (30) and (31) can he written in a more symmetrical foim if the 
solution is carried a step further. The equations for the set of con- 
stants that are used for a given configuration may he written, using 
the set o^s+l exan^le. 


-1 - Oi 7 oi + “3^03 + • • • “2s+1^0,2s+l 
^ ~ “l^U “3^13 ■*■•*• “ 2 s+l^l, 2 s+l 
0 = 0(^721 + ^3/23 + . . . «2 s + i 5'2,2 s +1 


(32) 


where the 7 's are the corresponding functions of and h. How 

consider as a matrix the coefficients of the ri^t-hand side of the 
equations whose constant terms are equal to zero, and let Ay 

A 5 , ... he equal to the determinants, respectively, which remain 
after the column which corresponds to the number 2 s + 1 of the con- 
stant A 2 S+I is removed. Then using Cramer’s rule, 


L 
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ai 


:^i 

^l701 - ^3^03 + • • • C-l)®-^s+l7o,2s+l 


^1/01 - A3703 


_^3 



(33) 


with, corresponding equations for the remaining values of Q2 b+.i. 

When the solution for each constant equation ( 33 ) is substituted 
into equation (30), it may he vrltten 


m/2 


dw 

35 - 


-phr(i - h^) IT - 2h^cos 29 - + 1 


g=i 


it(z2 - h^)(l 


m/2 


( 3 t) 


A?) U ^ - 2z^cos 20g + 1 
S =1 


The above equation must be multiplied by one of the following terms, the 
choice depending upon the symmetry of the desired slot configuration; 
for case I, or panels Intersecting both axes. 


it - 1] 

s=0 

2^ (-l)%3^lb2 ‘^®[^2(2s+1) ^ ^ 

s=0 


( 35 a) 


for case II, or slots intersecting both axes (note that in this case as 
well as in case IV, B has the same relation to p as A has to a 
in the demonstrated case). 


i 





^ (-l)%3,lb^""|>(2s4.l) . ^ 


( 35 b) 
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for case III, or 
the y-axls, 


and for case IV, 
the y-axis. 


a panel Intersecting the x-axls and a slot intersecting 
m+2 

(35c) 

(-l)®A2gb2 ^ i) 

a slot intersecting the x-axis and a panel intersecting 



m+2 

/ s I “ 2s+l 

1 2_ (-D^BggZ^ 

s=0 

m+2 

x: I-2S+1 

6=0 


(z**^ - l) 


- x) 


(35d) 


The constants A2s+l^ ®2s+l^ -^2s> ®2s determined from the 

following matrices in the manner discussed in the material following 
equation ( 32 ). The matrix for case I is 




cos 3© d© 



cos(| + l)© d© 
fl(0) 


(35e) 


r^m 

2 


-1 


cos © d© d© 


fi(o) 


2 " ^ 
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where 


m/2 

f2(e) = (l - 2b2cos 20 + h^) /cos^Qg - cos^e 

S=1 


Solution for symmetrically loaded wings .- Symmetrically loaded wings 
can he assumed, to he made up of vortex pairs of clrctilatlon -dr and dr. 
Since the circulation can he escpressed as 


r = r^f (x) (36) 

the strength of each pair ‘becomes 

dr=rof’(x)dx (37) 

This value of the circulation may he substituted Into eq.uatlon ( 3 ^) to 
give the contribution of each elemental vortex to the total flow. It Is 
then necessary to Integrate the equation over the entire span to obtain 
the flow, or 

where (<lw/dz)jj^ Is the complex velocity of the nonunlformly loaded wing, 
and (dw/dz)u. Is the ccorplex velocity (eq. (3^)) of the uniformly loaded 
wing with X substituted for h and a circulation of unity. 

Corrections for Interference dtie to lift .- The interference cnnglex- 
velocity field at the far-downstream position is determined by subtracting 
the con^jlex velocity of the free field frcm the ccsi 5 >lex velocity of the 
constrained field, or 


«= "as; ,(^2 . i2) 

A useful parameter which Indicates the effect of slotting a tunnel 
Is the ratio of the interference of that tunnel to the interference of 
a closed tunnel. This parameter, which is the same for the far-downstream 
position as at the model and which will be called the quality factor k, 
may be ejjpressed mathematically as follows; 
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itr 

^ * ,(z^ - 


f^S. + itr 

jr(z^ - 




where dw^/dz is the ccmplex velocity of the closed tunnel conf Iguraticn. 

Once the veilues of k: are detemilned as a function of the semi span b 

for a specific tiinnel, the interference may be calculated. It is shown 
in reference 9 that the interference on the lift may be expressed as an 
angle which may be added to the measured angle of attack to obtain the 
true or free-fli^t angle of attack. Reference 9 gives the correction 
angle in radians as 


= 


6SCl 

c 


ih-1) 


where 6 is a number which is determined from the geometry of the tunnel 
configuration and S is the wing area, C the tunnel area, and the 

lift coefficient. Also from reference 9, the increment which must be 
added to the measiired drag to obtain the correct drag is 

AOd = 6 ^ (^2) 


In the case of the slotted txumel. It is convenient to use the quality 
factor k and then express the two corrections as 


= k& ^ (i).3) 

SC ^ 

AQd = k5 — (44) 
o 

where the number 6 may be determined from the literature on lift inter- 
ference of closed tijnnelB. 

Solution and interference quality factors for tunnels with cross 
sections other than circular .- The solution of the wall- Interference 
problem for tvumels having cross sections other than circular may be 
obtained by following the method of paragraph 14.6, reference 10. This 
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method requires that a ftmction z = f(5) he found which will conformally 
transfoim the interior of the tunnel cross section in the z-plane into 
the interior of ttie unit circle |^j = 1 in the S-plane. It is also 
necessary that f'(^ does not vanish or hecome infinite in the unit 
circle | ^ | =1. The complex velocity of the txjnnel in the z-plane is 
then 


u - iv 


dw d^ 
d^ dz 


(45) 


In the function dw/d^ of equation (45), the 0n's which detenmine 
the slot edges of the ^ -plane are the transformed values of the slot edges 
in the z-plane and the values of Z in the ^-plane are also determined 
from the transformations of the points h in the z-plane. Once these 
values are determined for the transformed ttinnel in the ^ -plane, the 
velocity field dw/d^ in the ^ -plane may he confuted. The velocities 
in the z-plane may then he computed hy equation (45). 

The Interference and the quality factor may now he deduced frcm 
equation (45). Subtracting the free field frcm the cou^lex velocity 
given in equation (45) gives for the interference velocity: 


^ ^ + ihr 

dz dt; dz _ h^) 


(46) 


The interference for any tunnel can he determined from this equation; 
however, in many cases, the use of a quality factor may he more con- 
venient. The quality factor for this class of tunnels may he written 


dgdz ^(z2 - h^y 
^ ihr 

dz ^ jc^z2 _ h2) 


(47) 


where dw^/d^ represents the complex velocity of the transfoimatlon of 
the closed tunnel in the ^-plane. 

A simplification of eqtiatlon (47) may he made in case the transforma- 
tion f(5) may he approximated hy c^ for points near the center of the 
tvinnel. If h is sufficiently small, its value in the ^-plane may there- 
fore he represented hy h = cZ, wiiere Z represents the point which 
locates the transfoimatlon of the vortex in the ^-plane. Suhstitutlng 
these approximations into equation (47) gives for the quality factor: 
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Equation (48) shows that the quality fantor near the center of the 
original tunnel In the z-plane is approximately equal to the quality 
factor of the corresponding clrc\ilar tunnel In the ^-plane. Thus^ to 
obtain the interference, only the quality factor of the clrculeir tunnel 
needs to be ccm^juted, provided the interference of the original closed 
tunnel is known and the approximation z = cl is valid throoi^out the 
region in ■which the model is located. 


Applications of Theory to Various Tunnels 

Two slots located symmetrically across the x-axis (fig. 2(a)) .- 
Thls class of tunnels has the symmetry of slot location that is treated 
under case IV, so that its complex velocity field may be expressed by 
the product of equation (3^) and expression (35d) where s takes the 
values 0 and 1. Since the integrals across the two panels are identically 
zero, the matrix of the Integrals has no neaning and hence is not used. 
With these considerations the complex velocity may now be ■written 


(49) 

at(z2 - b^)(l - z2b^)/(z^ - 2z^cos 20i + l)B2(b4 - l) 

The ■wall interference is determined by subtracting the complex 
velocity for the free field from equation (49) ^ or 


ibr(l - b^)^b^ - 2b^cos + 1 - i) 



The quality factor k is determined by dividing equation (50) by the 
closed-tunnel interference which is 



J 
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or 


k = 


(l - z^)/b^ - 2b^cos Ze^ + 1. 

_ (l - (4^ - 2 z2cos 201 + 1 

(l - -b^)(l + z2) 

(l - z2b2) 


- 1 


(52) 


This f met ion may be written 


k = 


(l - a^)^b^ - 2b^cos 201 + 1 - (l - z2b2)^z^ - Sz^cos ZQj^ + 1 
(z^ - b^) ^z^ - 2z^cos 2 © 2 _ + 1 


(53) 


The effects of the slots on the Interference along the span of the model 
are obtained by substituting x for zj then. 


k = 


(l - x^)/b^ - 2b^cos 203 _ + 1 - (l - x^b^)^x^ - Ex^cos 202 . + 1 

- b^)^x^ - 2 x2cos Z6j^ + 1 


If X and b are sufficiently stnalT. so that the approximation 


(■ 


x^ - 2x^cos 202 . + 1 = 1- x^cos 20i + — sln^ 202 . 


is valid, then 


2 cos 20]_ - 2x^ - (x^ + b^)sin^ 202 ^ 
k = (55) 

2 - 2x^cos 202 _ + x^sin220i 

The qdallty factor is seen to be insensitive to valms of either x or 
b, so that for this tunnel with relatively small models (b = 0.25 or 
less) the effect of slotting the tunnel on the lift of any model may be 
obtained by multiplying the closed-tunnel interference for that model by 
the slotted- tiannel quality factor. 
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Pi stole si (ref. 7 ) shows that the Interference qtiality factor for 
this tunnel with a vortex doublet in the center is, in the notation of 
this paper, cos 202.. Equation ( 55 ) reduces to the same value \rtien x 
and h are set equal to zero. 

Two slots located symmetrically across the y-axis (fig. 2 (h)) .- 
This class of tunnels has the symmetry of slot location that is treated 
under case III (see eq. ( 35 c)). Thus, its ccnrplex-velocity field may he 
ejqDressed by the product of equation ( 3 ^) expression ( 35 c), so that 
it may be written 


dw 

dz 


-ibr(l - b^) |/b^ - 2b^cos 20^ + 1 

2Aqz 

- Ag(l + z^) 


jt(z2 - b^)(l - z 

2b2)j^ - 2 z2cos 201 

+ 1 

2 Aob2 - A2 (i + 

b“]l 


( 56 ) 


where Aq and A2 are determined frcan the matrix and are 


Aq 


cos 20 d0 


^-®1 (1 - Zb^coB 20 + b^)/cos% - cos202_ 


(57) 


'01 


d0 


J-02. (1 - Zb^cos 20 + b^))^cos20 - cos^0i 


( 58 ) 


If z and b are considered small the complex velocity of this 
field may be approximated by 


^^i ibr (^2^^ " ^* 0 ) 

1 + 1 

[1 . 

1 (Ag - 2 Aqz2] 

Icos 2©i 

<iz « (A2 - 2Aob2)(] 

L - z2b2) 


from which the quality factor becomes approximately 


(59) 


_ (Agb^ - 2 Aq) + (1 - b^) (Ag - 2 Aqz2) COS 202 . 


Ag - 2Aob" 


(60) 


uhlcli reduces for z = "b = 0 to 


k = - - — + COB 20 t (6i) 

As 

This quality factor is the same as the one presented in reference T, with the necesBary change 
of sign to adapt it to the notation of this paper. 

4s slots located symnetrlcally with respect to the x- and y-axes .- This general class of 
■h iiTiTiP lH Tirh thf» syirmetry of slot location of case If equation (^5a). if panels intersect hoth 
axes or of case II, equation (35b), if slots intersect both axes. Ihus, the coorplex velocity 
may he e3q>ressed by the product of equation (34) and expression (35a) or (35b), depending on the 
desired slot configuration. As an exan^ile, the coo^lex velocity for a tunnel containing four 
slots with panels intersecting both axes may be written 

-lljr(l - - a^eos 2Sg + l) + l) - + l)] 

^ !■ ; (62) 

k(z^ - b^){l - z2b2)fri l/z^ - SZ^COB 20g + lj|Aib^(b2 + l) - + 1^ 


where 


cos ^9 d8 

2 

(l - 2b^cos 20 + b^) /cos^0 - cos20j 


(63a) 
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and 


A 3 = 




- 0 -. 


cos 0 d0 


(63b) 


(i - 2t)^cos 20 + 13 ^) rr if 


cos^0 - cos 0 


g 


'1 g=l 

If the tunnel contains eight slots with panels intersecting hoth 


axes, 


dw 

dz 


-ihr(i - h^) XT - Sb^cos 20g + 1 

g=l 

5 ' 

it(z^ - h^)(l - XT ^z^ - -2 z2cos 20g + 1 

g=l 


A^z^^z^ + l) - A^z^^z^ + l) + As^z^® + l) 
A3^b^(b^ + l) - + l) + A5 (td 10 + l) 


( 64 ) 


where Ax^ A^^ and A5 ajre equal to the determinants remaining when 

the first, second, and third columns, respectively, are removed from the 
following matrix: 


P®1 cos 0 d0 

f®l cos 30 d0 

"'®l cos 50 d 0 

Le^ <■(«) J 

1 

CD 

H 

H> 

0^ 

f(e) 


i 


®3 cos 0 d0 

f(0) 


where 


'®3 cos 30 d0 

leg f(e) 


j; 


®3 cos 50 d0 

leg f(0) 


(65a) 


f (0) = (1 - 2b^cos 20 + h^) XT i/cos^0 - cos^^ 

g=l 


g 


(65b) 
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If the tunnel conta ins 12 slots with, panels Intersecting hoth. sjces^ 

6 

-ihr{i - h^) IT - 2h^cos 20g + 1 

X 


S =1 


dv 

dz 6 

- 2z^cos 20g + 1 

S =1 


A^^z^(z^ + l) - A^z^^z^ + l) + AgZ^^z^® + l) - A^^z^^ + ij 
Aih^(h^ + l) - A3h^(h^ + ^) + As'b^('b^° + 


( 66 ) 


where Aj^, A^, A5, and Aj are equal to the determinants remaining 

when the first, second, third, or fourth columns, respectively, are 
removed from the following matrix; 


r®l cos 0 d0 COB 30 d0 

;J-ei fCe) J-e^ f(«) 



COB 50 d0 

f(0) 


cos 76 d.0 

01 f(e) 


f^3 cos 0 d0 j 

cos 30 d0 

P ®3 COB 50 

d0 

J02 f(®) J 

02 ' 

Jsj f(s) 


f®5 cos 0 d0 j 

P ®5 cos 30 <10 

r®5 cos 50 

d0 

Jeif. f(0) J 

02^ f(0) 

J 04 f( 0 ) 




( 67 a) 



cos 70 <10 

f(0) 


where 


6 


* 

f ( 0 ) = (1 - Zb^cos 20 + h^) IT ^ 


fcos^0 


- COB<^0, 


g=l 


g 


( 67 b) 


If the slots are of equal width and evenly spaced, equations ( 62 ), ( 61 j-), 
and (66) can he simplified hy substituting the terms 


^1 - 2z™cos 201 + z^ 


and 


/l - 2b“cos 201 + 

for the multiple products in equations ( 6 g), (64)^ and ( 66 ), and 


LO 

o 


I 


9 m 9 m 

COS^ — 0 - C08^ ■g^ ©1 


for the multiple producte In equations (63a)> (63h), (65b), and { 67 b). With these suhstitutlons 

the limits of the Integrals in the second row of the matrix (65a) are changed to ~ - 0i for the 

4 

lower and ^ + 9l ^or the upper limit, and those of the second and third rows in the matrix ( 67 a) 
to ^ ■ ®1 ” ■ ®1 lower and ^ + ®1 f ^ upper limits. 

The Interference in these tumels is obtained hy coo^>utlng the Telocity and subtracting the 
free-field velocities, then dividing by I /2 since the coo^iuted fields Eire for a great distance 
downstream and (see r^, 8 ) the interference velocities at the model are one -half those down- 
stream. 

2(2r + 1 ) slots located symnetrlcally with respect to both aaes .- This general class of 
tunnels has the symmetry of slot location of either case III, escpresslon (35c), a panel inter- 
cepting the x-axis and a slot intercepting the y-axls, or case IV, expression (35d), which is 
the opposite of case III. Thus, the con^lex velocity may be expressed by the product of equa- 
tion ( 3 ^) aiTd expression (35c) or (35d), depeaidlng upon the desired configuration. As exan^les 
of this cIeibb and of the simplification due to evenly ^aced slots of equal width, Ihe coraplex- 
velocity fields for tunnels containing six evenly spaced slots of equal width may be written, 
for case III, panels on the x-axls. 


dv 

dz 


-ibr(l - b^l 

- 2b^cos 602 ^ + 1 

2Aqz 

- A 2 Z^(z^ + 1 ) + Ai(.(z® + 1 ) 

] 

jt(z2 _ b^)(l - _ Sz^cos 6d^ 

+ 1 

2Aob^ - Agb^Cb^ + l) + A4(b® 

j] 


( 68 ) 
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wiiere Aq, Ag, and. Aj^^ are evaliiated. from the following matrix hy 
the same methods used in equation (64): 


r«i de 

r®l cos 20 d0 

r®l cos 40 d0 

CD 

H 

CD 

f(0) J 

Le^ f(0) 


'1 + 01 


-- 01 


d0 

f(0) 




Jj-01 


COS 20 de 

f(e) 


^4 + 


01 


(69a) 


f-01 


cos 40 do 

" fCe ) 


where 


f(0) = (l - 2b^cos 20 + ■b^)\lcoB^39 - cos^30i 
For case IV, slots on the x-axis. 


(69b) 


dw 

dz 


where 


-ibr(: 

L - b^)03^2 _ 2h^cos 601 + 1 


[z'* - 1) - 24(28 - J.[ 

1 

*(22 - 

1^1 - - 2 z^cos 601 

+ 1 

B 2 b 2 (h^ - 1 ) - Bjj.(hS . 

■J 


(70) 


Bj, = 


— Q 

3 ®1 


sin 20 d0 


'01 (l - Bb^cos 20 + hM/cos2302_ - cos230 


(71a) 


and 


O 3-01 


Bg = 


sin 40 d0 


01 (1 - 2h^cos 20 + h^)^cos2302_ - cos^30 


(71b) 
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The Interferences are determined as in the previous sections. 

Equations (68) and (70) may he extended to 10, + 2 

slots in the same manner that equation (62) vas extended to inclvide 8 
and 12 slots. 


Single slot located symmetrically across the x~axis .- Since this 
tunnel c^talns an odd number of slots, sirnimation (I7 ) rather than (U) 
must he used in equation (l4), which, when expressed for a single slot, 
hecomes 

z y cos ^ ^ ^ ■6 + hjj sin 

dv _ n=0 ^ 

(z2 - h^)(l - z2h^)|/cos ■a - cos 

Equation (72) may he rewritten, with several modifications after sub- 
stituting z = e^'®, as 



dw 

dz 


a^z(l 


+ z) + a]_(l + z^) + i|^qz(1 - z) + hi(l - 
^ - h^)(l - z2b2)j^l - 2z cos + z2 



(73) 


The constants in this case must he evaluated from the circulation about 
the poles at h and -b and the nonexistence of sources at those points, 
as the potential condition is automatically satisfied because the, poten- 
tial is constEint across a single slot. The clrcvLLatlon and nonexistence 
of sources may be expressed about the pole b as 


r + i0 « 


2«i^b(l + b) + ai(l + b2) + l[bob(l - b) + bi(l - b^7|}‘ 
^ 2b(l - b^)^l - 2b cos 9 ^ + b^ 


(?ita) 


and atout the pole -b as 

2itl^o(-^)(l - "b) + ai(l + + l[bo(-b)(l + b) + bi(l - b3)]| 

,/2(-2b)(l - b^*-) )/l + 2b cos ©i + b2 


r + la = 


(T^t) 


If there are to be no soxirces at the poles, a must equal zero. It 
may be seen from inspection that ao and a^ are equal to zero and that 
bo and bi have solutions other than zero. When these solutions for 
bo arid hi are substituted into equation (73 )^ it becomes 



r 


*• y 


j 

r 




dv 

d2 


^(l + 1 d)/i - 2b GOB $]_ + + bCl - + 2b cos 0i + l)^z(l - z) + 

jl(l - b^} + 2b cofl 01 + b^ + ^1 + b3)^l - 2b cos + 1?](1 - .3) 


2b3c(z2 - b^)(i - zV)i/iirr 2z cos + z^ 


(75) 


The interference velocity Is determiiied as before by subtracting the free -field velocity 
-ibr/i(z2 - b^); hence, the interference . complex velocity is 

*"[b(l +b)^ 2b cos .0jL + b^ + b(l - b) /l + Zb cos + - .) + 


-ibr- 


dvj 

dT 


+s3)/r - 2b cos 0^ + b^ - (l - b^) ]/l + 2b cos 0^ + 1^5(1 . .3) . 
2b(l - " 2z COB + z^ 


2b«(z2 - b^)(l - zVj/lT' 2z COB 0]_ + 


( 76 ) 


The quality factor k is determined by dividing equation (76) by the closed-tunnel interference 

velocity, irtiich is ~ibr/ir(l ~ z^^) » After performing this operation and conibinlng a number of 
terms., the quality factor k may be written 


k != 3*. 


(1 + b)(l - z)(z + b)(l + zb)^l - 2b cos ©x + " 

(1 - b)(l - z)(z - b)(l - zb)^l + 2b cos + b^ - 
2b(l - z^b^)|/l - 2z cos 0x + 


2b(z^ - b2)/l" - 2z cos 03^ + z^ 


( 77 ) 


This factor may be checked by letting 0x = 0°, which represents a closed txmnel, for which k 

may be shown to be equal to +1; and if 0i = l8o°, idaich represents the open tunnel, then k can fjj 

be shown to be equal to -1. ^ 
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Equation (77) niay t>e consideratly reduced by using the following 
approximations, in which z is assumed to be small with respect to 1; 



2z cos 02_ + 


cos 02. + 


sin^02_ 


(78) 


/l + 2z cos ©2_ + = 1 + z cos ©2_ + z‘ 


3in20i 


(79) 


with similar eqiiations for the radicals containing b. After multiplica- 
tion and collection of terms, equation (77) reduces with the above 
approximations to 


sin^0n 


cos &i 


1 + b'^ 1 


k = 


b sin20. 


1 - Z cos 0]_ + z‘ 




( 80 ) 


If terms of the order z^ and higher are eliminated, eqiiation (80) may 
be approximated with 


k 


1 

2 


2 cos 02_ - cos%2. + 2 



- cos 



( 81 ) 


Eqioation (8l) shows that the interference of this timnel has an odd 
function con^jonent along the x-axis which will produce a variation in 
the effective angle of attack that will cause the model to have a 
rolling moment. Since extension of flow fields of this type (an odd 
number of slots symmetrically located with respect to the x-axis) to 
greater numbers of slots will not eliminate the odd powers of z, the 
rolling moments will continue to exist. It is therefore to be concluded 
that tunnels containing an odd mmiber of slots symmetrically located 
with respect to the x-axls Introduce an extraneous moment into the data. 

Single slot located symmetrically across the y-axis .- Since this 
t. unu pl Blan contninR n.ri niid mimb f»r nf fli nt.s, a ummfl. tinn (17) must be used 
to instore the proper slot symmetry. With this choice of functions, 
equation (15) may be written. 


dw 

dz 


Uq cos ^ + a^ cos -^ + bo sin ^ + b^ sin ^ 


(z2 - b^)(l - z^b^) 




sin d - sin 01 


(82) 
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Since z = the con 5 >lex velocity may he rewritten after suh- 

stitution and combining various terms as 

(l + i)[^z(l + z) + ai(l + z3) + holz(l - z) + hii(l - z3^ 

2 (z 2 _ h2)(l _ z2h2)/l - 21z sin Qj. “ 


(83) 


The constants &q, aj^, ho^ and h^^ must now he evaluated from 

the circulation about eanh of the poles and the condition of continuity. 
Ihe two conditions about each pole may he expressed as 


r + icT = 


2jti(l + i)jaQh(l + h) + ai(l + h3) + holh(l - h) + hil(l - h3^ 
l{-h(l - h^)\/l' - 2ih sin 02 _ - h^ 


(81m) 


about the pole h and 

2iti(l + i)[IaQh(l - h) + ai(l - h3) - hoih(l + h) + hil(l + b3j 
-4h(l - h2)l/l + 2ih sin Bi - h^ 

(81<.h) 


-r + ia = 


about the pole -b. The strength a must equal zero in order to satisfy 
the condition of nonexistence of sources. In order to simplify the 
rationalization of equation (84), ttie following substitutions are made: 



C + iD = ^ 

- b^ 

+ 21b sin 02 . 

and 






C - 

• ID = {l~ 

- b^ 

- 2ib sin 02 . 

Also let 






A-l 

= aQb(l + 

b) - 

f- a2.(l + b3) 


^2 

= -aob(l 

- t) 

+ a2.(l - b^) 



= bob(l - 

b) - 

K b2.(l - b3) 


®2 

= -bob(l 

+ t) 

+ bi(l + b3) 




and 


A = 


2br(l - b^)/l - 2b^cos 20i + b^ 


( 85 ) 
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The four simultaneous equations (84) may now te written, after ration- 
alizing the denominator. 


A = R.P.(1 - 
0 = I.P.(l - 
A = R.P.(1 - 
0 = I.P.Cl - 


i)(C + 1 D)(Ai + IBi) 
i)(C + 1 D)(Ai + IBJ 
i)(C - 1D)(A2 + IBg) 
i)(C - 1D)(A2 + IBg) 


which may be restated as 


A = Ai(C + D) + Bj^Cc - D) 

0 = -Ai(C - D) + Bi(C + D) 
A = A2(C - D) + B2(C + D) 

0 = -A2(C + D) + B2(C - D) 


If Ai and Bi are determined from the first two and A2 and B2 
from the last two. It may be seen that 

-gT = Ai = aob(l + b) + ai(l + b^) 

2(c2 + 

~ = A2 = -aob(l - b) + ai(l - b 3 ) 

2(c2 + D^; 


fio “ = Bl = + ^1(1 - 1 =^) 


21 + D 


7V '~~ 2) " = -bo^Cl + T3) + tl{l + b3) 

2(c2 + D^) 
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These equations may he solved for sq, a^, hQ^ and hj_, giving 


&Q - 


a.n — 


hn = 


A 


(c + 

D)(l - h3) 

- (C - 

D)(l + 

h3) 

+ 



lt-h(l - 

■hM 



A 


(C + 

D)h(l - h) 

+ (C - 

D)h(l 

+ h) 

C2 + 

d2 


lrt)(l - 

■ hM 




A 


+ D^(l - h3) - fc 

Di(l 

+ t3) 


+ E 

,2 

i4h(l 

-hM 



A 


(c + 

D)h(l - h) 

+ (c - 

D)h(l 

+ ^>) _ 


r (89) 


C2 + 


4h(l - hM 

The set of equations ( 89 ) may he simplified to 


an = -hn = 


» 2a(d - Cb8) 

4(c2 + D2)h(l - h^) 


r (90) 


ai = hi = 


2Ah(C - Ifti) 

l4-(c2 + ri2)b(l - h^) 


Letting aQ = -hQ a n d ai = hi, the complex velocity equation ( 83 ) may 
he rewritten 


dw 

dz 


2[z(ao + s^z^) + l(ao“'^ + 

S(z 2 - 1 , 2 ) (1 . a 2 b 2 )l/r - 2 iz sin 01 - z 2 


(91) 
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SulDstituting for eq and a^ gives 


dv _ 
dz 


2A j 


[p - Cb3 + b(C - Db)z3 

+ 1 

[b(C - Db) + (d - Cb3)z^ 

L 

k(c2 + D2)b(l - b^)(z2 - b2) 

(i- 

z2b2))/l - 21z sin 02 ^ - z 

2 


(92) 


Now It may be shown that 

A = - ^(l - b^)(c2 + D^) 


(93) 


so that 

M = - £ 

dz Jt 
where 


;|(d - Cb3) + ^(c . iJb)z^ + l|b(C - Db) + (d - Cb3)z^ 
(z2 - b2)(l - z2b2)/r - 2iz sin 0i - z2 


(9lt) 


C + ID = /l - b^ + 2ib sin 

The interference velocity is obtained subtracting -ibr/rt(z2 _ b2), 
SO tlaat 

l)(c - Db) + (d - Cb3)z2 _ 

|1z[d - Cb3 + b(C - Db)z^ 


r~ (K 


IbF 


dw^ 

«(z2 - b2) 




b^l - z2b2)^l - 2iz sin 0]_ - z2 


- 1 


(95) 


Now divide by the closed- tunnel interference to obtain the qu a l i ty factor, 
so that 

b(C - Db) + (d - Cb3)z2 - iz |d - Cb3 + b(C - Db)z2] - 


k 


•s 

-ibr 

b(l - z2b2)l| 

t - 2iz sin ^ J 

Jt(z2 - b2) b(l - z2b2)i 

^1 - 21z sin 0]^ - z2 


-ibr 


(l - b2 

)(l + z2) 

1 

jt(z^ - t 

.2) 

1 - 

z2b2 


( 96 ) 
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or 


}s. = 


b(C - Db) + (d - Cb3)z2 _ 1 z|d - Cb3 + b(C - Db)z^] ^ 
b(l - - 2iz sin ©i - z2 


^(z^ - b^)/l - 2iz sin ~ ^ 


(97) 


The Interference quality factor at the center of the tunnel may he 
examined hy letting z = Oj then 


k = 


C - Dh - 1 

Ih^ 


(98) 


Frcm the evaluation of the constants, equation (9^), 


when 


\^1 - 2h^cos 2©i + h^ cos t 

( 99 ) 

1 - Zh^cos 2©]_ + h^ sin 

( 100 ) 

- , 2h sin ©T 


+ . I tan-l ^ 

( 101 ) 


When h is small, equation ( 98 ) may he simplified hy first approxi- 
mating C D. These approximations are, for smal T values of h. 


2 

C = 1 - ~^cos 201 + sin20l) 


( 102 ) 


D = h sin ©1 ( 103 ) 


When the values for C and D given in equations (102) and ( 103 ) are 
substituted into equation ( 98 ), it hecomes 


cos^©! 


k = sin ©1 + 


2 


(10k) 
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This equation also checks with the results given in reference 7^ 
with the usueO. change in sign to conform with the notation of this paper. 


Lift interference of rectangular tunnels .- The function which 
transforms a rectangular tunnel into a circular tunnel is given in 
paragraph l^i-.S, reference 10, as 


5 = 


cn is. 
2 


(105) 


Xz Xz* 

where sn — , dn — , and cn — are Jacobian elliptic fimctions of z 
2 2 2 

and x/z is defined by 


X ^ K ^ ^ 
2 a h 


(106) 


where K and iK' are the qviarter periods and a and h are the 
breadth and, height of the tunnel. (See pars. l4.7 and 14.8 of ref. 10 
for fiarther information concerning these functions.) 

If a and h are. given, then K, K', and m, the squared modxilus 
(see ref. 10, par.' l4.8), are uniquely determined. Once the constants K, 
K' , and m are determined, the slot location and the half -span length Z 
may be coiBputed. First, consider slots which are located on the top and 
bottom of the tunnel (see fig. lb). In reference 10, paragraph 14.8, it 

is shown that the top of the tunnel may be expressed by z = x + ~ or 

Xz = ?tx + iK' . It is also shown that the transformation (l05) may also 
be expressed 


^2 ^ ( 107 ) 

1 + cn Xz 

so that the equation for the top of the tunnel becomes, in the C plane, 

f2 ^ 1 - cn(Xx: + iK' ) 

1 + cn(Xx + iK' ) 

From reference 10, paragraph 14.8, cn(Xx '+ IK') is equal to -m?-/2ds iXx:, 
which may be called i|i. Now, on the circle. 
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f 2 ^ g210 ^ J- - Ipt 

^ 1 + in 

and if tbe real and imaginary portions are separated. 


cos 20 = 


1 . - -tif 
1 + 


m - (ds Xx)^ 
m + (ds 


(109) 


(no) 


or 


sin 20 


2u _ 2ds Xx 1 
, 2 ~ .f— ^ . 6lb ~Xx 

1 + M- Vm 1 + — 


( 111 ) 


The ©ri's of the circular tunnel in the ^ plane are determined from 
the above equations by the locations of the slot edges Xn in the z-plane. 

If the slots are on the side of the tunnel, a similar analysis 
stairting with z = ^ + iy or \z = K + iXy will show that 


sin 20H 


2\lm sd iXy 
1 + m(sd ily)^ 


( 112 ) 


As before, equation (H2) may be used to determi ne the 0n's when the 
slot edges are functions of y. 


The value of I is given in reference 10 as 


1 - cn Xb 
1 + cn Xb 


(113) 


Knowing I and the 0n's, the con^jlex velocity of the tunnel in the 
§ plane may now be coniputed, and the interference determined from 
equation (46) . 

The ceilculation of the interference for rectangular tunnels may be 
simplified by applying the conclusion following equation (48), which 


42 


NACA m L53A26 


states that the quality factor of the transformed tunnel is approximately 
equal to the quality factor of the original tunnel, provided the trans- 
formation function can he approximated with z - c^. Equation (l05) is 
shown in reference 10, paragraph 14.8, to he approximately equal to 
^ = Xz/2; hence the conclusion following equation (48) is valid for 
rectangular tunnels* Thus, it is necessary only to determine the quality 
factor for the transformed tunnel and applying it to the correction for 
the fully closed rectangular tunnel (see ref. 10, par. 14.8, for this 
correction) to obtain the correction for the corresponding slotted rectan- 
gular tunnel. 

Circular tunnels having general slot configurations .- Equation (l6) 
is used to obtain the complex velocity for a configuration which contains 
a general slot distribution. Equation (l6) may be written 


dw 

dz 


m ^ 

pz(a^^ cos -d + b^^ sin <^ 3 ) ^ (o^ + 1 ( 311 )^ 0 ^^ 2 

n ?=0 



(11^) 


where p has the same definition as given for equation (25) . If z Is 
suhstituted for ■9 by uBlng the relation z = e^^, equation (ll4) becomes 


pz 


dw 

dz 


(*1 ^ (“n ife) (l f^)“ 

2m ^ 

n,k - cot ^ 

g=lN 21-1 2 


(115) 


Equation (ll5) may be reduced to 


dw 

dz 


p|a2^(l + z^) + ibi(l - z2^ 



(ctQ + ipn)(z - + L)“i“ 




z2b2 


^ 

) rr v^(z + 1) - (z - i)cot ^ 

g=i ’ ^ 


(116) 


0 


« 


Ik 


The constants oc> “l; • • • °n# 3o^ Pl^ • • • Pn determined in the 

same manner as vas used for the symmetrical slot configurations; that is, the circulation about 
both poles is equal to F, the source strength a is zero, and the potential in each slot is 
equal to zero. The conditions on the circulation and source strength may be determined from the 
line Integral about the poles b and -b. For the pole at b. 


r + ia 


2itip[a3_(l + b^) + lbi(l - b^l] ^ (“n + iPu)(l - b)”^^(l + b)^l’^(-l)“"“ 


2oi 


I 

4b (l - b^) riwi(b + 1 ) - (b - l)cot 

g^l 


2 


(117a) 


and for the pole at -b^ 


F + Iq B 


2itlp ax(l + b^) + lbx(l - b^J (on + lPn)(l + ^r"’^(l - 'b)’^^(-l) 

2m I 

4b(l - b^) n \/l(l - b) + (1 + b)cot ^ 
g=l * ^ 


m-n 


(UTb) 


The condition on the potential in each slot may be determined from integrating the velocity over 
each panel, or 


0 = 


^ a 

^a^ cos 0 + bx sin 0) > (ax + iPx)cot“ - 
n^ 

[l - 2b^cos 20 + - cot ^ 


d0 


(U-Tc) 


£ 


panel 
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The constants a^; hi; oq, ai, . . . OnJ and Pq> • • • 3n 

are determined from the solution of the set of simultaneous equations (117). 
Once the constants are determined, the- complex velocity may he determined, 
Sind the Interference may then he computed In the same manner as was used 
for the various symmetrical cases. 


RESULTS MD DISCUSSION 


Interference quality factors for several, circular tunnels- with 
symmetrically located, evenly spaced slots .- The quality factors for 
circular tvmnels containing various symmetrically located, evenly speiced 
slots of equal length and wings of very small span (h = O) are given In 
figure 3* curves for tunnels that have panels across the x-axls and 

contain 4, 8, or 12 slots are calculated from tlie formula 


k = - 


Jo 






e de 


cos2 “ e - cos2 




e de 


I I e - cos2 


(118) 


where 6 Is defined In figure 1 and m Is the number of slots. Eqioa- 
tlon (ll8) may he derived hy considering the value of eqioatlon (34), 
expressed for symmetrically located, evenly spaced slots of equal width, 
when h and z are equal to zero. It Is also the negative of the 
relation given In reference 7 for the same tunnel configurations, the 
minus sign helng used to conform with the notation of this paper. 

An analysis of figure 3 shows that, for all configurations except 
the one with two slots across the x-axls, only a sinal 1 percentage of the 
t unn el wall must he opened In order to obtain no Interference at the 
center of a Rmen model, and that the amotint of opening required rapidly 
becomes smaller when the number of slots Is Increased. It Is also noted 
that the change In quality factor Is very rapid at the null-interference 
condition, thus making an accurate estimate of the Interference dlfflcialt. 

The Ingjortance of slot location In determining the qtiallty factor 
Is Indicated hy the large divergence of the two curves for tunnels con- 
taining two slots. It can he shown, though, hy comparing the qimllty- 
f actor functions for the conditions of slots across the x- and y-axes 


NACA RM L53A26 


k5 


and panels across tbe x- and y-axes, that as the nvnnher of slots Increases, 
the quality factors will approach each other. Since the quality factors 
approach each other, it may he expected that the quality-factor curves 
for tunnels containing 8 or 12 evenly spaced slots of equal length will 
be approximately cozrrect for any slot location. Thus, the quality factor 
should be approximately correct even thou^ the model is rolled in the 
tunnel. 

Variation of quality factor with span of model .- The effect of model 
span on the quality factor k at the center of several slotted tiamels 
is shown in figure i^. Examination of the ciirves shows that the quality 
factor is fairly constant over an appreciable range of spans, in that 
it does not vary more than 10 percent for spans of 0,5 to 0.6 of the 
tunnel diameter or width. The quality factor for the tunnel containing 
two slots is seen to vary more, in that the quality fantor changes 0.10 
for a change in the span from zero to 0.46 of the tunnel diameter. The 
quality factor of the single slot, however, is seen to remain fairly 
constant regardless of span. 

The spanwise variations of the quality factor for a model with a 
span of 0.75 of the tunnel diameter is shown in figure 5. This figure 
shows that the quality factor for the tunnels containing 8 or 12 slots 
is about 0.3 larger at the tips than it is at the center, whereas for 
the tunnel containing two slots it decreases about 0.6. These changes 
indicate that the spanwise change in k may depend tqion the boundary 
condition at the Intersection of the wall and x-axisj that la, if a 
panel intersects the x-axis, k will become larger, and if a slot inter- 
sects the x-axis, k will become smaller. Siich a tendency would be 
difficult to prove, however, so that an analysis of the variation in k 
for any particular tunnel would req^lire a ccmqrutatlon of k along the 
span for that model. These variations in k also indicate that the 
lift corrections for a model with a span as great as 0.75 of the tunnel 
diameter can be roughly approximated by using the value of k at the 
center of the tionnelj however, a more accurate correction would involve 
the use of an average value for k as well as an average of the load 
over the span. 

A T^Buai& of the observations made from figures 4 and 5 indicates 
that, for these specific slot configurations at least, the quality 
factors for tunnels containing 8 or 12 slots, provided they are symmet- 
rically located, equally spaced, and of equal length, are more nearly 
constant throughout a greater portion of the central region of the tunnel 
than are the quality factors in tunnels containing two slots. The fact 
that the quality factor is moderately constant even at spans of 0.75 
permits making the correction for lift Interference by conqputlng the 
Interference of the closed tunnel having the same cross section, and 
then multiplying that interference by the quality factor for the center 
of the slotted tunnel. 
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Quality factors In •b\miiels containing a single slot .- The variation 
of the quality factor with percent opening in a tunnel containing a 
single slot symmetrically located with respect to the x-axis is shown 
in equation ( 80 ) and in figure 6 to be somewhat different from, those 
previously studied in that k is greater on the half -span of the wing 
which points toward the panel and is smaller on the side which points 
toward the slot. This variation of k means that a spanwlse variation 
in angle of attack exists which would cause the model to roll. Thus, a 
tunnel of this type, that is, a single* slot symmetrically located with 
respect to the x-axls, has interferences which are more difficult to 
correct for because of the introduction of an unnecessary rolling moment 
into the measured data. 

This rolling moment does not exist, however, if the slot is symmet- 
rically located with respect to the y-axls. It is shown in equaticm (97)^ 
thou^, that the induced velocities will not be nomnal to the span but 
will have a small cang>onent along the span. Such a component should not 
affect the total forces seriously, since it would not affect the total 
induced velocity. However, the effects of this cco^onent should perhaps 
be considered in the treat^nt of load distributions along the span. 

Thus, the pecularltles of the interference effects of tunnels containing 
a single slot symmetrically located with respect to either axis indicate 
that its effects will be more difficiilt to correct. 

If the equations for the single-slot case are extended to a larger 
odd number of evenly spaced, equal-width slots, it may be expected that 
both the rolling moment and the cross flows will become smaller because 
the walls create a more uniform interference field at the model. 

Effects of ccatpressibillty on the corrections .- It is shown in 
reference 11 that the lift due to cotpressible flow can be corrected 
in exactly the same manner as thou^ the flow were incompressible. Thus, 
neither the k nor the 5 of equations (k3) and (Wt) is a function of 
compressibility. Since the arguments used in reference 11 are based on 
subsonic linearized conipresslble flow, it may be expected that correc- 
tions can be made for Mach nunibers tap to the critical or sli^tly higher. 

Applicability of the theory .- Several of the differences between 
the idealized and the actuEil paroblem are -Uiose due to viscosity, which 
causes a mixing region in the neighborhood of ~Uae slot mther than the 
assumed constant pressure surface which divides the hi^-veloclty tunnel 
air and the stagnant tank air. Since the mixing region involves various 
ccoplicated phenomena such as turbulence, velocity gradients, separation 
at the outer slot edges, and differences in boundary conditions in dif- 
ferent slots depending upon the dlrectlcmi of flow throu^ the slots, it 
is very likely that its effects on the lift interference will have to be 
determined from analysis of experimental data (for example, ref. 12). 
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Another effect of viscosity •which heccmes Ingjortant if the slots 
are narrow nnt^ deep is the friction on the air as it flows throu^ the 
slots. Since the friction reduces the amovint of flow throu^ the slots, 
they ■will effectively heccane narrower, so that the quality factor may 
be expected to become larger. 

The possibility also exists that, in an actual tunnel, the pressures 
may not be equal in all the slots and the slot pressures may not be equal 
to the pressiire in "the tank. It would be expected that if the difference 
between the pressxire in the slot and that in •fahe tank causes more flow 
throu^ the slot than occurs in the ideal tunnel, hhe slot will effec- 
tively be wider, and if the difference decreases the flow, the slot will 
effectively be narrower . Thus, the actixal quality factor ■will depend 
upon the effective slot ■wld'ths as determined by the pressxire increments. 

Other differences occvir because a practical t\mnel cannot be con- 
structed like the ideal tunnel. These differences involve finite slot 
lengths, 'variable-'width slots, and lips on the slot edges. In con- 
sidering the effects of finite slot lengths, it seems reasonable to assume 
that those effects should be no more serious than the effects of a finite- 
length open tunnel. In reference 13^ it is shown that the theoretical 
lift corrections for an infinitely long open tunnel are adequa-te, pro- 
vided the model is located a distance of at least one -half the tunnel 
helgjht frcan the entrance and exit regions. Therefore, if the slot con- 
figuration is STich that the -width is constant over a section whose length 
is at least equal to tunnel height, the tlieoretlcal corrections shoiold 
be adequate e'ven thou^ the slot width may -vary considerably outside that 
region. The effect of lips on the slots may be shown qualitatively by 
comparing the press\ire gradients of the two types of slot configurations, 
that is, ■with and without lips. Since the lips confine the flow, the 
pressure gradient in that configuration is less steep through the slot 
than if there were no lips. Therefore, the velocity may be expected to 
be lower, so that "the effective -width of -the slots -will be reduced, 
thereby increasing the quall'ty factor. 


COlKLlIDIISrG REMARKS . 


An analysis of the equations -which represent the interference on 
the tralldng-vortex system of a uniformly loaded wing, due to -wind- 
tunnel walls with mixed open and closed boundaries, has shown that: 

1. Slot openings of the order of 7 percent for four e-venly spaced 
slots of equal length, and less for larger numbers of slots, are required 
to reduce the interference on a lifting model to zero. The zero- 
interference quality factor is critical with respect to the percent of 
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slot opening as a small change in slot opening will cause an appreciable 
change in the quality factor when its value is near zero. 

2. The tunnels which contain two symmetrically located slots showed 
guite different values of the Interference for different slot locations. 
The differences in interference became sms, Her as the number of slots 
increased. 

3. In the tunnels examined, a region existed about the center of 
the tunnel in which the ratio of the slotted-tunnel interference to the 
closed-t\mnel interference was fairly constant, so that in order to 
obtain the corrections for the effects of a slotted tunnel it is neces- 
sary only to multiply the closed-tunnel interference by a constant. 

4. The region in which the ratio of the slotted-tunnel interference 
to the closed-tvuanel interference is reasonably uniform was fo\md to be 
larger for the tunnels containing 8 or 12 slots than .for those containing 
2 slots. 

5. An examination of tunnels containing a single slot showed that 
they produced a rolling moment or a cross flow on the model. These 
phenomena may be expected for tunnels containing a larger odd number of 
evenly spaced, equal-width slots. These phenomena should, however, 
decrease as the number of slots increases. These phenomena may also be 
expected for any case in which the slots are asymmetrically arranged 
with respect to the model axes. 

6. An analysis of the effects of coirpresslbllity on the ratio of the 
interference of a slotted tunnel to the interference of a closed tunnel 
shows that the ratio is relatively unaffected by compressibility through- 
out the subsonic region. 


Leuigley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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(b) Square tunnel . 


Figure 1.- The coordinate systan used for the investigation of the lift 
interference due to slotted tunnels. 
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Slots Percent opening 
in tunnel wall 
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Ratio of distance from center of tunnel 
to tunnel radius, x/R 


Figure 5 . - Variation of the quality factor along the span in several 
tunnels containing equally spaced, symmetrically located slots of 
equal length. Span of model is 0.75 of tunnel diameter. 
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Point at which interference Slot 

is calculated configuration 

0 ^ 



Percent of opening in tunnel wall 


Figure 6.- Quality factors for a tunnel containing a single slot and a 
model whose span equals one-half the tunnel diameter. 

UNCLASSIFIEu 


NACA-Langley - 4-?-53 - 92S 



tJNCLASSiFlED 


tf^JCLASSinFO 


